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We summarize our recent results on the resummation of hard-scattering co- 
. efficient functions and on-shell form factors in massless perturbative QCD. 

The threshold resummation has been extended to the fourth logarithmic 
order for deep-inelastic scattering, Drell-Yan lepton pair production and 
Higgs production via gluon-gluon fusion. The leading six infrared pole 
^ il terms have been derived to all orders in the strong coupling constant for 

D . the photon-quark-quark and the (heavy-top) Higgs-gluon-gluon form fac- 

' tors. These results have many implications, most notably they lead to a 

new best estimate for the Higgs production cross section at the LHC. 

X 

5h ' 1. Introduction 

Coefficient functions, or partonic cross sections, form the backbone of 
perturbative QCD. These quantities are calculable as a power series in the 
strong coupling constant as, but exhibit large logarithmic corrections close 
to threshold. The all-order resummation of the dominant soft-gluon contri- 
butions takes the form of an exponentiation in Mellin-A^ space [1-4], where 
the moments are defined with respect to the appropriate scaling variable, 
like Bjorken-x in deep-inelastic scattering (DIS) and x = M^+^_ ^/s for the 
Drell-Yan (DY) process and Higgs production via gluon-gluon fusion. 

The purpose of the exponentiation is (at least) two-fold. On the one 
hand, it can directly lead to improved phenomenological predictions close 
to exceptional kinematic points, for instance to an improved stability un- 
der scale variations. On the other hand, it can be viewed as a generating 
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functional of fixed-order perturbation theory close to the partonic thresh- 
olds. Hence progress in the soft-gluon resummation also facilitates improved 
fixed-order predictions which, depending on the specific observable, can be 
relevant even very far from the hadronic threshold. 

In this contribution we discuss recent results for the threshold resumma- 
tion up to the fourth logarithmic (N^LL) order [5,6], and briefly illustrate 
their implications. We also summarize our recent results [7, 8] for the on- 
shell quark and gluon form factors and their exponentiation [9-12], which 
were instrumental in extending the soft-gluon resummation to N^LL accu- 
racy for lepton-pair and Higgs boson production. Moreover the form-factor 
results are interesting also in a wider context, e.g., they provide another link 
to recent calculations performed in AA = 4 Super- Yang-Mills theory [13]. 

2. General structure of the threshold resummation 

As mentioned in the introduction, the coefficient functions for inclusive 
DIS, Drell-Yan lepton-pair production and Higgs boson production expo- 
nentiate after transformation to Mellin A^-space [1,2], 

= (1 + as 501 + «s' 502 + • • •) • exp (G^) + 0{N-^ In" N) . (1) 

Here collects the A^-independent contributions at A;-th order in the strong 
coupling constant Us- The resummation exponent contains terms of the 
form In^ A^ to all orders in and takes the form 

= \nN-g^{\) + 52(A) + 53(A) + 54(A) + ... (2) 

with A = /3o Os In A^. The functions 5^ represents the contributions of the 
fe-th logarithmic (N'^~^LL) order. All our relations refer to the MS scheme. 

The exponential in Eq. is build up from universal radiative factors 
A p and Jp due to radiation collinear to the initial- and final-state partons, 
and a process-dependent contribution A"^* from large-angle soft gluons. For 
example, the resummation exponents for the processes considered here read 

G^iS = InAq + InJq + InAg^g , 

CfoY.H} = 2 1nA{q_g} + lnA|DY,H} • (3) 

Ap, the so-called jet function Jp and A'°* are given by certain integrals over 
functions of the running coupling, A^, Bp and D. Specifically, the functional 
dependences are Ap(^p), Jp{Ap,Bp) and A™'^(Z)). The functions Ap, Bp 
and D, in turn, are defined in terms of power expansions in ag, for which 
we generally employ the convention 

k=i 
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The extent to which these functions are known sets the accuracy to which 
the threshold logarithms can be resummed. It is worth noting that the 
function -Dj^jg is found to vanish to all orders [14,15], hence Aj^jg = 1. 

The explicit expressions for the functions (7j(A) in Eq. Q are obtained by 
performing the above-mentioned integrations, for instance using properties 
of harmonic sums and algorithms for the evaluation of nested sums [16-19]. 
Specifically, and g'4 have been determined in Refs. [20,21] and [5], to which 
the reader is referred for details. While the leading-log (LL) function 
depends only on Ai, the N'^-^LL functions gi._^_i include all parameters up to 
Afc+i) -Sfc and Z?fc. We now turn to the present status of their determination. 

3. The known resummation coefficients 

The functions Ap are given by the leading large- (or large-x) coeffi- 
cients of the diagonal splitting functions for the parton evolution. 



Ppp(as) = Ap(as) (1 - + Pp\a,) 6(1 - x) + O (ln(l - x)) 



(5) 



which in turn are identical to the anomalous dimension of a Wilson line with 
a cusp [22]. The known expansion coefficients for the quark case read [23,24] 
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for rij effectively massless quark flavours. Here Cp and Ca are the usual 
colour factors (Ci;' = 4/3, = 3 in QCD), and Riemann's zeta function is 
denoted by Cn- The gluonic coefficients are related to Eqs. © by [22,25] 



A, 



Ca/ Cf ^q,i . 



(7) 



It is worthwhile to note that the C,2 terms in 3 have been confirmed by 
the recent 7\A=4 Super- Yang-Mills (SYM) calculation of Ref. [13]. 

The perturbative expansion of the functions Ap{as) is very benign. In 
fact, already A3 has a very small effect on the resummed coefficient functions 
[20,21]. Therefore it is sufficient to estimate the presently unknown fourth- 
order coefficients A^ entering g^ by their [1/1] Fade approximants. 



Aq,4 ^ 7849, 4313, 1553 for n. = 3 , 4 , 5 



(8) 
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to which we assign a conservative 50% uncertainty in numerical apphcations. 
Eqs. © and © lead to the numerical four-flavour expansion 

A^{as,nf = 4) ^ 0.4244as (1 + 0.6381as + 0.5100a,2 + 0.4 [i/ij + . . .) • (9) 

We now turn to the coefficients Bp entering the jet functions Jp. These 
quantities can be determined by comparing the ag-expansion of Eqs. 
and with the results of fixed-order calculations of the DIS coefficient 
functions, which we have recently extended to the third order in [26]: 
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The result for Bq^i is, of course, well-known [1,2], and Bq^2 has been derived 
by us before in Ref. [27] where we explicitly established also D 



DIS 



0. 



For the extraction of Bq^^ [5], on the other hand, we rely on the all-order 
proofs [14, 15] of -Ddis = mentioned above. 

The numerical expansion of Bq in QCD is far less stable than Eq. @, 



5q(as,nj = 4) 
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Note, however, that the large third-order contribution to Bq actually sta- 
bilizes the expansion of shown in Figd for Bq^s = and iV = 40, for 
example, the N^LL term would be about as large as the previous order. 

The coefficients Bg i for the gluonic jet function Jg are, for instance, rel- 
evant in direct-photon production which is dominated by the qq — > gj and 
qg —f q7 subprocesses close to threshold, see Ref. [28]. These coefficients 
can be obtained in the same manner as Eqs. (jlOj) . but from DIS by exchange 
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Fig. 1. Left: successive approximations for the resummation exponent ^ of inclu- 
sive DIS. Right: minimal-prescription [3] convolutions with a typical input shape. 

of a scalar (j) with a pointlike coupling to gluons, like the Higgs boson in 
limit of a heavy top quark. We have derived the corresponding coefficient 
function C(^,dis ^^P to the third order in the course of calculating the lower 
row of the flavour-singlet splitting function matrix [25] . Comparison of these 
results to the expansion of Eq. yields -Bg^i and the previously unknown 
quantities Bg^2 and -Bg,3- The analytic results can be found in Ref. [5]. Here 
we confine ourselves to the numerical expansion in four-flavour QCD, 

Sg(as,n^ = 4) ^ -0.6631 as (1 - 0.7651 - 2.696 + . . .) , (12) 

which shows a third-order enhancement similar to that in Eq. (|llj) . 

Finally we address the process-dependent coefficients Di due to the large- 
angle emission of soft gluons. Up to now, the two-loop coefficient functions 
for proton-proton processes are known only for the Drell-Yan cross section 
and Higgs boson production in the heavy-top approximation [29-32]. The 

corresponding coefficients D^^'^^ have been extracted from these results in 
Refs. [20,21]. Even for these processes, the three-loop coefficient functions 
have not been calculated so far. It is possible, however, to derive their third- 
order coefficients from mass- factorization constraints [6] , using our recent 
results for the pole terms of the three-loop quark and gluon form factors [7,8] 
and the third-order splitting functions [24,25]. Postponing the discussion 
of this derivation to section 5, the results for DY case read 
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(13) 



The corresponding coefficients for Higgs boson production via gluon-gluon 
fusion are found to be related to these results by a simple colour-factor 
substitution, 



Df = Ca/Cf dI 



DY 



(14) 



which is in complete analogy to Eq. Q. It worth pointing out that both the 
cusp anomalous dimensions and the coefficients and exhibit a 
maximally non-abelian colour structure, as anticipated for in Ref. [22]. 

The numerical expansion of in four-flavour QCD is given by 



Z)^^(as,nj = 4) ^ 2.3211 as (0 + + 2.675 < + . . .) 



(15) 



The ratio of the third- and second-order coefficients is very similar to that 
for the jet function in Eq. (|lip. underlining the numerical relevance of D^. 



4. On-shell form factors and their exponentiation 

The form factors of quarks and gluons are gauge invariant (but infrared 
divergent) parts of the perturbative corrections to inclusive hard scatter- 
ing processes. They summarize the QCD corrections to the qqX and ggX 
vertices with a colour-neutral particle X of either space-like or time-like 
momentum q. These quantities are also key ingredients in the infrared fac- 
torization of general higher-order amplitudes [33,34]. 

The relevant amplitude for the space-like ^*qq case is 

= iCq (u7^u) J'^q(as,Q^) , (16) 

where Cq represents the quark charge and = —q^ the virtuality of the 
photon. The gauge-invariant scalar function ^q is the space-like quark form 
factor which can be calculated order by order in the strong coupling in 
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dimensional regularization with Z) = 4 — 2e. The corresponding Hgg vertex 
defining is an effective interaction in the limit of a heavy top quark, 

£eff = --CnHGI^^G"-'^" , (17) 

where G^^^ denotes the gluon field strength tensor, and the prefactor Ch 
includes all QCD corrections, known to N^LO [35], to the top-quark loop. 

The well-known exponentiation of the form factors is achieved by 
solving the evolution equations [9-11] 



Q^TT^lnJ^ as,— ,e = - K^a^.e) + - G [ —,as,e\ (18) 



based on a factorization of the form factor T into two functions K and 
(jr. The latter are subject to renormalization group equations [9] which are 
both governed by the same anomalous dimension A-p of Eqs. © and 
because, obviously, the sum of G and K in Eq. (|18j) is a renormalization- 
group invariant. We follow the decomposition of Refs. [11,36], where the 
function K \s a, pure counter-term collecting the infrared 1/e poles, while 
the infrared- finite function G includes all dependence on the scale . 

The resummed form factor is given as a double integral with the bound- 
ary condition ^(as,0,e) = 1 [11]. After both integrations are performed, 
\n.J- exhibits double logarithms of Q'^ j ^ and double poles in e. The rela- 
tion H18() can be then used for a finite-order expansion and matching of the 
predictions to the results of explicit higher-order calculations. The resulting 
expressions for the bare expansion coefficients T\ in terms of the quantities 
Ay and the (still e-dependent) Og-expansion coefficients Gi of G(Cf' j [i^ = 1) 
in Eq. (|18l) are sketched below (see Ref. [7] for the complete formulae): 

= 8^^' + - ^o) + h^Gl + ...-A2)-^G2 

^3 = -4^A? + ... + ^(9GiG2 + ...-4^3)-^G3 

•^4 = ^vlf + ... + ^(3G2 + 8GiG3 + ...-3^4)-^G4 . (19) 

We have extracted all three-loop pole terms of the quark and gluon form 
factors J-c^ and J-g from the calculation of the third-order coefficient functions 
for DIS by the exchange of a photon (coupling to quarks) and a scalar 
(coupling to gluons) [26], already mentioned above in the discussion of the 
jet function Jp. The details will be reviewed in the next section. 
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Similar to the two- loop analysis of Ref. [12], we write the coefficients Gp as 
Gp,! = 2(^P^i- 5pg/3o) + /i' + eGp^i , 
Gp,2 = 2 (Pp% - 2,5pg/?i) + fi + /?oG;;i(e = 0) + eG;^2 , 
Gp,3 = 2 (Pp% - 3,5pg/?2) + fi + PiG;',i{e = 0) 

+ f3o[G^,2{^ = 0) - PoG^,i{^ = 0)] +eG;',3 (20) 

with F = [F — F(e = 0)]. The quantities Pp have been defined in 
Eq. (ISJ above, and the terms with Jpg are due to the renormalization of the 
operator G^j^yG^'^ in Eq. ((T7|l . The crucial point of the decomposition (PU)) 
is that the functions turn out to be universal and, like the Ap in Eqs. © 
and (O maximally non-Abelian with (at least up to the third order) 

ft = Ca/Cf ft . (21) 

The explicit results for the quark case read 

ft = 0, ft = 2C7^|-/?oC2-^n^ + Q(^-14C3) 

_ ^ ^2 / 136781 12650 1316 . , 352 . , 176 
/s -CfCa[^-^ ^C2-^C3 + — C2 + — C2C3 + I92C5 

11842 2828 _ 728 . 96 ^2 1711 



+ CaCfu, + ^ C2 + ^ C3 - y C2^j + C^n, 

304 ^ 32 .\ ^ 2^ 2080 40 ^ 11-^, 
+ 4C2 + -C3 + y Clj +C.n/ (-^ - ^C2 + ^C3j . (22) 

Note that ft has been obtained already in Ref. [12], and that the coefficients 
of the highest ("-function weights, C2C3 Cs &t three loops, agree with the 
results inferred from the recent AA = 4 SYM calculation in Ref. [13]. 

Going back to Eq. ()19() . it is worth noting that the leading term of G3 
in Eq. (|20|) . together with corresponding coefficients of Gi and G2 to higher 
powers in e (see Refs. [7, 8] for the explicit results) fix the six highest poles 
of the form factors at four loops and, in fact, at all higher orders. Moreover, 
taking into account that the numerical effect of A4 in Eq. Q is small, our 
present results are sufficient for deriving the infrared finite absolute ratio 
|^p(g^)/^p(— (7^)P of the time-like and space-like form factors up to the 
fourth order in as- The corresponding numerical results for = 4, 5 read 

qq'j* : 1 + 2.094 Os + 5.613 + 15.70 + (48.63 ± 0.43) , 
ggH : 1 + 4.712 Os + 13.69 + 25.94 + (36.65 ± 0.35) , (23) 

where the the uncertainty of the last terms is due that of the fourth-order 
cusp anomalous dimensions ^p,4, as estimated below Eq. (jSJ in section 2. 
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5. Partonic cross section and their infrared pole structure 

In this section, we finally discuss the extraction of the form factors from 
our calculation of the coefficient functions for inclusive DIS and the related 
derivation of all soft-enhanced third-order terms for the Drell-Yan process 
and Higgs production, and thus of given already in Eqs. (|13j) and (|14|) . 
from these form- factor results and mass-factorization constraints [6]. 

The starting points for the first step are the explicit results for the 
bare (unrenormalized and unfactorized) partonic structure functions for 
7*q — > qX and (p*g ^ gX in the limit x — > 1 [26]. At each order keep- 
ing only the singular pieces proportional to (5(1 — x) and the +-distributions 



V, 



(1-x) 



/ = 1, . . . 2n - 1 , (24) 



these results are compared to the general structure of the n-th order con- 
tribution F^ in terms of the /-loop form factors J-i and the corresponding 
real-emission parts 5/ , 

^0 = '5(1 -x) 

F^ = 2Ti6{l-x) + Si 

F^ = (2.F2 + .Fi^) 5(1 - x) + 2JPi5i + 52 

F^ = {2J^3 + 2TiJ^2) S{1 - x) + (2J^2 + J^i) Si + 2J^iS2 + S3 . (25) 

In DIS the x-dependence of the real emission factors S^ is of the form 
Skifk,e)i with the D-dimensional -|— distributions fk^e defined by 

A,,(x) = 6[(l-x)-i-^^]+ = -l5il-x) + Y^t^eV,. (26) 

1=0 

The dimensionally regularized (with D = 4 — 2e) bare structure functions 
F^ in Eq. (|25() exhibit poles in e up to e"'^", with a structure completely 
determined by mass factorization. On the other hand, the individual real 
and virtual contributions Tk and Sk in Eq. (|25|) contain poles up to order 
e~2A:^ which cancel due to the Kinoshita-Lee-Nauenberg theorem [37,38]. 

The determination of the form factor now proceeds as follows. Once the 
combinations of lower-order quantities in Eq. (|25() have been subtracted from 
F^, the n-loop form factor can simply be extracted by the substitution 

Vo ^ Ls{l-x)-Y.t^v,, (27) 

ne ^ I ! 

1=1 
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which exploits the particular analytical dependence oiSk on x, i.e., Eq. (|^ . 
As 6{1 — x) enters with a factor 1/e, this extraction loses one power in e. 
Hence from the third-order calculation to order e'', as performed for the 
coefficient function, we can only extract all pole terms of .F3 in this manner. 

The second step, the determination of the +-distribution contributions 
to coefficient functions for lepton-pair and and Higgs boson production, 
proceeds along similar lines, see Ref. [39] for an early two-loop application 
to the Drell-Yan process. In analogy to Eq. (|25|) . the soft limit of the bare 
partonic cross sections for qq — > 7* — > and gg H reads 

< = S{l-x) 

W}' = 2ReJ^i5(l -x) +5i 

= (2Re^2 + |^i!^)'^(l-x) + 2Re^i5i+52 

= (2Re^3 + 2|.Fi^2|)5(l-x) + (2Re.F2 + |^i|^)5i 

+ 2Re.Fi52 + 53 , (28) 

where, of course, J-" now denotes the time-like quark or gluon form factor, 
known by analytic continuation from = —Q^ < to > 0. The real- 
emission contributions Sk depend on the scaling variable x = M^, h/^- 
this case, the dependence of Sk on x is of the form Sk{f2k,e)i i-e., 

00 

Sk = hk,e E 2A:Sfc,,e' . (29) 

l=-2k 

With the known time-like form factors, the expansion coefficients ^ of 
the soft function Sk can be derived recursively as far as they are subject 
to the KLN cancellations and the mass-factorization structure relating the 
remaining poles to the splitting functions Q. Employing the results of 
Refs. [7,8] and [24,25], the third-order terms S3 _g ... S3^_i can be obtained. 
Due to Eq. 1)26(1 this is sufficient to derive all +-distribution contributions to 
the third-order coefficient functions, in particular also the coefficient of T>q 
from which D^^^'^^ in Eqs. (tTTni and (fill) , can be determined by matching. 
An important application on these new results is presented in Fig. |21 

The connection between mass-factorization and resummation leads to a 
simple relation between the coefficients and /p in Eqs. (|21|) and (|22() . 

^^{DY,H} ^ _2J2 + 2/3oSi_o 

^^{DY,H} ^ _2/3 + 2/3i5i,o-4/3o'si,i + 4/3o[ 



*"2,0 ^ Ci ^{F,A} 



,(30) 



which has also been derived by extending the threshold resummation to the 
A^-independent contributions [40,41], see also Ref. [42]. In our approach, 
the I terms can be traced back to the ag-renormalization of Eqs. ()28() . 
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Fig. 2. The perturbative expansion of the total cross section for Higgs production 
at the LHC. Left: dependence on the Higgs mass A/h. Right: renormahzation-scale 
(in-)stabihty for Mh = 120 GeV. See Ref. [6] for a detailed discussion. 



6. Summary 

Building on our third-order computation of the splitting functions [24,25] 
and the coefficient functions for inclusive DIS [26], we have derived new 
three- loop and all-order results for the threshold resummation [5,6], the 
on-shell quark and gluon form factors [7,8], and the coefficient functions 
for lepton-pair and Higgs boson production at proton colliders [6]. These 
results have important implications within and beyond perturbative QCD. 
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